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Hopping expansion as a tool for handling dual variables in lattice models 

Christof Gattringer a * 

a Massachusetts Institute of Technology, Center for Theoretical Physics 
77 Massachusetts Avenue, Cambridge MA 02139 USA 

The hopping expansion of 8- vertex models in their Grassmann representation is studied. We use the functional 
similarity of the Ising model in this expansion with the hopping expansion of 2-D Wilson fermions to show that 
the lattice fermions are equivalent to the Self-avoiding Loop Model at bending rigidity l/y/2. 



1. Grasssmann representation for 8- vertex 
models and the hopping expansion 

Many lattice models have representations in 
terms of dual variables. E.g. the partition func- 
tion of the 2-dimensional Ising model can be writ- 
ten as a sum over all possible Peierls contours 
with a weight given by e~ 2 ^ to the power of 
the length of the contours. A generalized model 
which includes many dual models of lattice sys- 
tems is the 8-vertex model In the 8- vertex 
model the dual variables (loops) are decomposed 
into their elements and each of those elements 
(vertices) is assigned a weight. The vertices can 
be viewed as 8 different tiles and are shown in 
Fig. [fl. The figure also gives the weights u>i in the 
notation of [fLl . 
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Figure 1. Vertices and their weights in the 
vertex model. 



The loops are obtained as tilings of the plane 
such that there is never an open end for the fat 
line. The partition function of the 8-vertex model 
is then given by 

z *> = e n (^) n< . 

loops 2 — 1 

where u>i is the weight of the i-th vertex, and ni 
gives the abundance of this vertex in the loop 
configuration. 

It is known that the 8-vertex model has 

a representation as an integral over Grassmann 
variables 



8r 



= J [dry] e 
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(2) 
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Here to each site x of the dual lattice a vector of 
4 Grassmann variables 

T)(x) = (r 1+1 (x),n- 1 (x),'n + 2{x) 7 r]_2{x)) T , (3) 

was assigned. The action S is a sum of propaga- 
tor-, corner- and monomer-terms (S = S p + S c + 
S m ), with 

s p = E t a V+i{x)V-i{x + i) 

X 

+b ry +2 (a;)7 ? _ 2 (a; + 2)] , 
S c = E t c y+i(x)v~2(x) + d 77 +2 (x)?7_i(a:) 

X 

+e T}-2(x)l]-i(x) + f f] + 2{x)f] +1 (x)} , 
S m = E [ 9 V-l(x)V+l(x) + h T]-2(x)T} + 2(x)].(i) 

X 

When expressing the weights uii through the co- 
efficients a, b, ... h as 

u>i = —cd — ef + gh , L02 = —ab , L03 = bg , 
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co>4 = ah , ^5 = cVab . loq — dVab 
LU7 = eVab , ujg — fVab , 



(5) 



it can be shown by explicit expansion that the 
Grassmann integral (0) reproduces the partition 
function ([!]). Since (^) can be diagonalized by 
Fourier transform, finding a Grassmann represen- 
tation (0), (|J) corresponds to solving the model. 
This is reflected by the fact, that the choice (||) 
obeys the free fermion condition CO1LU2 + uj^lu^ = 
which is a sufficient condition for an 
analytic solution of the 8- vertex model pj . 

Anti-symmctrization turns the action into a 
quadratic form S = \ J2 x . y v( x ) T K( x , y) 1 l(y) 
with kernel 



±2 
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M 



E p » 



(G) 
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where M, P^ have lattice indices x, y and indices 
i,j acting on the Grassmann 4- vectors These 
matrices are given by 



M{x,y) 



and 

(P+l(x,y))i 



Jx,y 
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(7) 



^-1 = -Pit 



P-2 = -Pi 



(8) 



Besides the termwise expansion which gives the 
partition function (Q), the Grassmann integral (^) 
can also be evaluated as a Pfaffian giving 



Z Hv = PfK = VdetK 



(9) 



In the second equality we used the fact that K 
is anti-symmetric and thus its Pfaffian reduces to 
the root of a determinant. The matrix M has 
determinant (gh — ef — cd) 2 which after inserting 
(H) and using the free fermion condition reduces 
to det M = uj 2 . Thus for 8-vertex models with 
lui =/= 0, the matrix M can be inverted and one 
can transform to 



Z 8v = VdctAiVdetll + ff] = 



\ n— 1 



(10) 



where we used the well known formula for the 
determinant as the exponential of the trace of the 
logarithm and defined the hopping matrix H as 



H 



±2 

E 



(11) 



The series in the exponent converges for values of 
the parameters where || H ||< 1. Since the ma- 
trices P M contain the hopping factors 8 x+il>y only 
terms that correspond to closed loops can con- 
tribute to (|l(i|). These closed loops are however 
related to the loops used in the original definition 
([!]) of the model and thus the exponent in is 
in principle the free energy of the model expressed 
in terms of loops. The only remaining problem is 
to explicitly solve the traces over powers of the 
hopping matrix H in (|l0|). For the general case 
this will be discussed elsewhere [|) . 

The Ising model can be obtained from the 8- 
vertex model by setting 



1 , w 2 



-2/3 



J =3, ...8. (12) 



For this case the hopping expansion (n0[) gives 



exp 



E 



-20\l\ 



(13) 



Here I runs over all closed, simply connected loops 
on the lattice, n(l) denotes the number of self- 
intersections, 1(1) gives the number of times the 
loop I iterates its entire trace and \l\ is the length 
of the loop. 

2. Equivalence of 2-D Wilson fermions 
with the Self-avoiding Loop Model at 
bending rigidity l/\/2 

In this section we give an example how repre- 
sentations of the type ( |l0| ) can be used. We will 
show that 2-dimcnsional Wilson fermions on the 
lattice are equivalent to the Self-avoiding Loop 
Model at bending rigidity 1 / For an alter- 
native derivation of this result see Scharnhorst's 
paper §. 
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The action for 2-dimensional Wilson fermions 
(§) is given by 



S = 4)Mil) 
with 



±2 

/i=±i 



M 



(14) 



(15) 



where r ±M = ~[1 =p /x = 1,2. Here cri, cr 2 
denote the Pauli matrices and the hopping pa- 
rameter k is related to the bare fermion mass m 
via K = (m+2) _1 . Since the action for the lattice 
fermions is a bilinear form the partition function 
is not a Pfaffian but a determinant 



Z w (k) = /[#][#] e 



detM 



/ J^r n 

det[l - kQ] = exp - ^ — TrQ r 



(16) 



Again the contributions to the traces in the ex- 
ponent correspond to closed loops, and the re- 
maining problem is to compute the traces of the 
matrices T M ordered along the links of the loops. 
This problem was solved in || by realizing that 
the Pauli matrices give rise to a representation of 
rotations on the lattice. The resulting hopping 
expansion is 



Z w = exp ( 2 k' 



1 



(_1)«(0 
~lW~ VV2 



c(i) x 



(17) 



This result is very similar to Formula ( |13| ) for the 
Ising model, with the only difference of an extra 
weight l/y/2 for each corner of a loop I (here c(l) 
denotes the number of corners in a loop I). 

When expanding the exponential functions in 
( fl3"| ) or ( |l7j ) one obtains contributions from prod- 
ucts of loops I appearing in the exponent, such 
that in these products some links are multiply 
occupied. In addition already single loops I can 
retrace some of their links such that they are 
multiply occupied. For the Ising model we know 
that all the contributions with multiply occupied 
links cancel since the Peierls contour representa- 
tion discussed in the beginning allows only for 



simply occupied links (see || for a detailed dis- 
cussion). This cancellation mechanism is inde- 
pendent of corner weights since two loop config- 
urations which differ in their number of corners 
necessarily have different occupancy numbers at 
some links and thus only contributions with equal 
numbers of corners can cancel each other. Hence, 
using the similarity between (|lj]) and (^) we 
can conclude that also in an expansion of the ex- 
ponential in (p"7|), terms with multiply occupied 
links cancel. The remaining terms inherit the cor- 
ner factors l/y2, and the self-intersection factor 
(— l) n W in ( p~7| ) leads to a cancellation of terms 
with self-intersections (i.e. the corresponding ver- 
tex model has weight u>2 — and thus the loops 
are self-avoiding). Thus we end up with (see |^] 
for the detailed proof) 

V Z w (k) = Z S al m (k, , (18) 

with the Self-avoiding Loop model defined as 
Z salm (z,v) = 5> |7| ?? C(7) ■ (19) 
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Here 7 runs over all self-avoiding contours, c(-y) 
gives the number of corners, r\ is the bending 
rigidity, z is the monomer weight and \ j\ denotes 
the length of the contour. 
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